II 62 50608 



CASE FILE 
COPY 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



REPORT No. 608 



STRESS ANALYSIS OF BEAMS WITH SHEAR 
DEFORMATION OF THE FLANGES 




1937 




REPRODUCED BY 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

U.S. DEPARTMENT OF COMMERCE 
SPRINGFIELD, VA. 22161 




AERONAUTIC SYMBOLS 

1. FUNDAMENTAL AND DERIVED UNITS 



Length. 
Time__. 
Force__ 

Power_. 
Speed. _ 



Symbol 



I 
t 

F 



P 
V 



Metric 



Unit 



meter 

second 

weight of 1 kilogram. 



horsepower (metric), 
fkilometers per hour. 
\meters per second. _ 



Abbrevia- 
tion 



m 
s 

kg 



k.p.h. 
m.p.s. 



English 



Unit 



foot (or mile) 

second (or hour)__ 
weight of 1 pound 

horsepower 

miles per hour 

feet per second 



Abbrevia- 
tion 



ft. (or mi.) 
sec. (or nr.) 
lb. 



hp. 

m.p.h. 
f.p.s. 



2. GENERAL SYMBOLS 



g, 

m, 

i, 

Mi 

8, 

S w , 

G, 

b, 

c, 

b 2 

S' 

V, 

L, 

d, 

Do, 

D u 

D„ 

C, 

R, 



Weight =mg 

Standard acceleration of gravity = 9. 80665 
m/s 2 or 32.1740 ft./sec. 2 

A1 W 

Mass= — 
g 

Moment of inertia=ra& 2 . (Indicate axis of 

radius of gyration k by proper subscript.) 
Coefficient of viscosity 



v, Kinematic viscosity 

p, Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-nr 4 -s 2 

15° C. and 760 mm; or 0.002378 lb.-ft " 4 sec. 2 
Specific weight of "standard" air, 1.2255 kg/m 3 

0.07651 lb./cu. ft. 



at 



or 



3. AERODYNAMIC SYMBOLS 



Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio 

True air speed 



Dynamic pressure=2P^ 2 

Lift, absolute coefficient C L =^ 



D 



Drag, absolute coefficient C D =-^ 
Profile drag, absolute coefficient Cd 0 =^ 



Induced drag, absolute coefficient C^—^g 



Parasite drag, absolute coefficient @d p = ^§ 

Q 

Cross-wind force, absolute coefficient C , c == ^§ 
Resultant force 



0, 

VI 
p — > 



Ola* 



Angle of setting of wings (relative to thrust 
line) 

Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 

Reynolds Number, where I is a linear dimension 
(e.g., for a model airfoil 3 in. chord, 100 
m.p.h. normal pressure at 15° C, the cor- 
responding number is 234,000; or for a model 
of 10 cm chord, 40 m.p.s., the corresponding 
number is 274,000) 

Center-of-pressure coefficient (ratio of distance 
of c.p. from leading edge to chord length) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 
Flight-path angle 
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STRESS ANALYSIS OF BEAMS WITH SHEAR DEFORMATION OF THE FLANGES 

By Paul Kuhn 



SUMMARY 

The fundamental action of shear deformation of the 
flanges is discussed on the basis of simplifying assumptions. 
The theory is developed to the point of giving analytical 
solutions for simple cases of beams and of skin-stringer 
panels under axial load. Strain-gage tests on a tension 
panel and on a beam corresponding to these simple cases 
are described and the results are compared with analytical 
results. For wifig beams , an approximate method of 
applying the theory is given. As an alternative, the 
construction of a mechanical analyzer is advocated. 

INTRODUCTION 

The so-called "semiinonocoque" type of const met ion, 
which has been favored by aircraft designers for some 
time, presents serious difficulties in stress analysis. 
Static tests have proved that the bending action of such 
a structure is not always described with sufficient 
accuracy by the standard engineering formulas based 
on the assumption that plane cross sections remain 
plane. It will be necessary, therefore, to devise new 
working theories for the action of semimonocoque beams 
under bending loads. 

In order to arrive at reasonably rapid methods of 
stress analysis, it is necessary to make rather sweeping 
assumptions. It is obvious that the range of applica- 
bility of any such method is limited. The present 
paper concerns itself with beams typical in general 
form of one class of beams used in airplane construction, 
that is, with fairly shallow, wide beams, having flat covers, 
symmetrical about the center line, with two shear webs and 
with bulkheads that offer no appreciable resistance to 
deformation out of their planes. 

Briefly, the action of such a beam underloads applied 
at the shear webs is as follows: The transverse shear is 
taken up by the shear webs. The flanges attached to 
these shear webs furnish part of the longitudinal stresses 
required to balance the external bending moment. 
The strains set up by these stresses induce shear stresses 
in the skin which, in turn, cause longitudinal stresses in 
(he intermediate stringers attached to the -kin until 
sufficient longitudinal stresses exist at any section to 
balance the external bending moment. 

If the skin between stringers did not deform under 
the action of the shear stresses, the standard beam 
formulas would apply. The thin sheet, however, has 



very little shear stiffness and suffers large deformations 
under load. As a result, the first intermediate stringer 
next to a shear web carries a smaller stress than the 
flange of the shear web, the next intermediate stringer 
carries less stress than the first one, and so on to the 
center stringer, winch carries the smallest stress. This 
phenomenon of the interdependence between stringer 
stresses and shear deformations forms the subject of 
the present paper. 

Apparently Dr. Younger was the first person in this 
country to give serious attention to this subject. In 
reference 1 he gives a formula for the efficiency of a box 
beam with walls of uniform thickness, which may be 
considered as the limiting case of very many extremely 
small stringers. Nothing more on the subject was 
published until two experimental studies appeared in 
193G. Reference 2, dealing with the case of a skin- 
stringer panel in edge compression, includes n t heoret ical 
solution for a particular case. Reference 3 deals with a 
box beam in pure bending, a problem identical with the 
one treated in reference 2. In both studies the stringer 
stresses experimentally obtained were used to compute 
efficiency factors for the shear stiffness of the sheet. 

The most important practical problem is the inverse 
of the problem dealt with in references 2 and 3 ; namely, 
given the shear stiffness, to calculate the stringer 
stresses. The problem is difficult and complex. In 
order to arrive at any solution, it has been necessary to 
use a very much simplified concept of the action of the 
structure, as suggested in references 1 and 2. On the 
basis of this simplified concept, the analytical solutions 
for a few very simple cases of axially loaded panels and 
of beams are derived in this paper. For other cases, 
it will be shown that a trial-and-error method of solution 
is feasible. 

The analytical solutions as w r ell as the trial-and-error 
method apply only to very elementary cases, namely, 
to three-stringer panels under axial load and to beams 
with a single longitudinal stringer attached at the 
center line of the cover sheet. It has been considered 
worth while to devote considerable space to the dis- 
cussion of these elementary cases for the following 
reasons: 

1 . The study of these simple cases greatly facilitates 
the understanding of the fundamental principles. (It 
is very strongly urged that anyone desiring to use the 

1 



2 
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proposed method of analysis work, by the trial-and- 
error method, at least one example each of a panel under 
axial load and of a beam.) 

2. The simple cases afford a very convenient way of 
experimentally checking the validity of the assumptions 
made. Strain-gage tests made for this purpose on a 
tension panel and on a beam are described in this paper. 

3. The solutions obtained for beams with a single 
longitudinal can be used as checks on the degree of 
approximation attainable with the " constant-stress 
method" proposed later for analyzing actual wing beams. 

An additional reason for the lengthy discussion will 
only be mentioned in passing. Under certain condi- 
tions, a beam with a single longitudinal stringer may 
give useful approximations of the stresses in a beam 
with many stringers. Such a simplified substitute beam 
makes it possible to obtain some rough ideas on the 
influence of bulkheads, an influence that was neglected 
in the present discussion. 

Two methods are proposed for winglike structures. 
One method is the construction of a mechanical ana- 

f 




(a) Axially loaded panel subjected to shear deformation. 

(b) Mechanical model. 

Figure 1.— Three-stringer panel. 

lyzer permitting a solution that is "exact" within the 
assumptions made. The other method is based on the 
assumption that the structure is so dimensioned as to 
approach the ideal design of constant flange stress along 
the span. For this ideal case, the analytical solution 
can be obtained. 'Die actual case will have deviations 
from the ideal case, which are termed "faults." These 
faults are minimized as much as possible by applying 
corrections, and the stresses caused by the corrections 
are superposed on the stresses of the ideal case. 

SYMMETRICAL THREE-STRINGER PANEL UNDER AXIAL 

LOAD 

FUNDAMENTAL CONSIDERATIONS 

The simplest possible structure in which shear 
deformation must be taken into account is shown in 
figure 1 (a). Two stringers, 4 1 and A' , of equal section, 
are connected to an intermediate stringer B by means 
of a thin sheet C. The upper edge of this sheet is 
reinforced by bars D. The stringers and the sheet are 
:ii i ached to a foundation l . 



The important phases of the elastic action of this 
structure may be visualized with the help of the 
mechanical model sketched in figure 1 (b). This 
model represents one-half the structure, which is per- 
missible because the structure is symmetrical. Helical 
springs represent the stringers A and B and their 
elastic resistance to longitudinal deformation. Coil 
springs represent the elastic resistance of the sheet to 
shear deformation. It is assumed that the stringers 
carry only longitudinal stresses and that the sheet 
carries only shear stresses. For the mechanical model 
it is assumed that guides prevent any deflection of the 
springs other than that for which they are designed. 

The stresses resulting from the load P are shown 
qualitatively in figure 2. At the top of stringer A the 

LP 




Figure 2.— Notation for axially loaded panels. 

stress is a A =P/A A , at the top of stringer B it is 0-5 = 0. 
The shear stresses r acting on the sheet gradually take 
the load out of stringer A and transfer it to stringer 
B. If the panel has sufficient length and if the sheet 
has sufficient shear stiffness, the stresses a A and a n will 
be very nearly equal at the root. 

EQUATIONS OF THE PROBLEM 

The equations governing the problem under the 
simplifying assumptions can be very easily set up. 
Figure 3 shows a strip of length dx cut from the panel 

F A +dF A F s + dFB 



T 



dS c 



\ dS c 



B 



dx 

I 



r 



Figure 3.— Element of panel. 



and separated into its component parts. 'The equation 
of equilibrium gives 



dF A = dSn — — dF B 



(1) 



(Sec list of symbols, appendix A.) 

It should be noted that these equations are written 
for the structure as shown in figures 1 (b), 2, and 3, 
w hich is one-half the original structure in figure 1 (a), 
so that A B is one-half the area of stringer B as shown in 
figure 1 (a). The sign convention used throughout this 
paper is that tensile forces and stresses are positive and 
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that shear forces and stresses in the sheet are positive 
when caused by positive stresses in the loaded stringer 
A (or in the flange F in the case of beams). 

The elastic deformation of the structure is shown in 
figure 4. Two corresponding points 1 and 2 are dis- 




^^^^^^^^^^^ 

Figure 4.— Elastic deformation of panel. 

placed to new positions 1' and 2' '. The total displace- 
ments are given by 

Ua = Jo ii dx and Ub== jo 

The shear strain is given by 



and since 



u A —u B 



7=77 



where G e is the effective shear modulus, these relations 
may be combined into 



rbE f* , p , 

-75—== I (Tidx— I ov/.r 
b e JO JO 



The last equation may be written 

(/r = ^(a A — a B )dx 



(2) 



Equations (1) and (2) may be combined into a differ- 
ential equation (see appendix B) which, together with 
the boundary conditions, defines the problem com- 
pletely. If there are more stringers, a system of 
simultaneous differential equations results. 

SOLUTION OF THE EQUATIONS 

For the fundamental case of a symmetrical three- 
stringer panel of constant cross section, the analytical 
solutions are given in appendix B for two cases: The 
panel attached to a rigid foundation and loaded at the 
free end, and the panel free in space strained by displac- 
ing the ends of the stringers a known amount. Com- 
bining the two solutions makes it possible to calculate 
loaded panels attached to an elastic-ally yielding 
fo 1 nidation. 

For the analysis of three-stringer panels in which the 
stringer areas and the shear stiffness of the sheet vary 
along the axis, a brial-and-error method has been found 
feasible. 



The recommended procedure for the trial-and-error 
method is as follows: 

Divide the length L of the specimen into a suitable 
number of bays. Tabulate the average values of t, A A , 
and A B for each bay. 

Assume values for the increment of shear AS C in 
each bay. According to equation (1) 

AF A = -AF B =AS C 

With the assumed values of AF A and AF B and the 
known values F A =P and F B =0 at the end of the panel, 
calculate for all stations along the length of the panel 
the forces in the stringers and then the stresses in the 
stringers. From these values calculate the shear stresses 
and the shear forces in the sheet. The method of tabu- 
lation is shown in table I. In this example, the values 
of A A , A By and t are constant and need not be tabulated. 

The calculated values of AS C will not, in general, agree 
with the originally assumed values. Change the assumed 
values and repeat the entire process until a satisfactory 
agreement is reached between the assumed values of 
AS C and the calculated ones. 

In the choice of the first set of values for AS C , the 
analyst must be guided by previous experience. The 
only condition known at the outset is 



S c < 



A A +A B 

because this is the maximum possible force that would 
be transmitted to stringer B only if the shear deforma- 
tion were reduced to zero. 

The most difficult step, and the one upon which the 
success of the method hinges, is to compare the cal- 
culated AS C curve with the assumed one and, on the 
basis of this comparison, to derive a new curve modified 
in such a way that the repetition of the entire calculation 
will yield a calculated AS C curve that agrees with the 
assumed one. No general rule can be given concerning 
the method beyond stating that decreasing the assumed 
AS C values at any point will raise the calculated ones and 
vice versa. Some practice is necessary to develop the 
skill required for this step. Five trials should be suffi- 
cient, in general, to obtain an agreement to 1 or 2 per- 
cent for five or six bays unless the variations of areas are 
extreme. 

It should be emphasized that the method is a trial- 
and-error one and not a method of successive approxi- 
mation, i. e., the calculated AS C curve cannot be used 
as the assumed curve for the next cycle. 

EFFECTIVE SHEAR STIFFNESS AND EFFECTIVE STRINGER AREAS 

Two quantities must be determined before an analysis 
can be started — the effective shear stiffnesses and the 
effective stringer areas. 

The shear stiffness of a flat sheet is equal to the shear 
modulus G of the material. If the sheet buckles into a 
diagonal-tension field and the edge members are rigid, 

the shear stiffness is the theoretical shear stiffness of a 

diagonal-tension field G € t=%G (for duralumin or steel 
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Sec tion A- A 
Figure 5.— Test panel. 



The condition of a pure diagonal-tension field is not 
reached, however, until the buckling shear stress has 

been considerably ex- 
ceeded. Consequently, 
values intermediate be- 
tween G and %G will 
occur at stresses not 
too greatly in excess of 
the buckling stress (i. e., 
3 to 5 times), provided 
that the edge members 
are sufficiently stiff. If 
the edge members are 
not sufficiently stiff or 
well braced to take the 
t ra nsverse component of 
the diagonal tension and 
particularly if the sheet 
ea fries edge compression 
in addition to shear, the 
shear stiffness may drop 
to very low values. 
Values as low as G e = 
0. Iff have been reported 
(reference 3); although 
the numerical accuracy 
of this particular anal- 
ysis has been questioned, it serves at least as a useful indi- 
cation of what may be expected, remembering that this 
test was stopped long before 
reaching the ultimate load. 
Quantitati v e information o n 
this subject is scarce. Fortu- 
nately, as will be shown later, 
the shear stiffness need not be 
very accurately known to obtain 
reasonable accuracy in the 
stringer stresses. 

It is clear that the sheet will 
not only act as a shear member in 
accordance with the theory but 
will also assist in carrying longi- 
tudinal stresses. The following 
assumptions have been used: 

1. For a sheet carrying ten- 
sion in addition to shear, it was 
assumed that the sheet is full} 7 
effective in tension; i. e., the 
sheet up to a line hallway be- 
tween the stringers is added to 
the stringer proper when com- 
puting the cross-sectional area 
of the stringer. This assump- 
tion is obviously somewhat un- 
safe and should be modified when 
the stringer stresses are high. 

2. For a sheet carrying com- 
pression iu addition to the shear, 



von Kar man's 
in the form 



formula for effective width was used 



where w is the effective width (on one side of the 
stringer) and <r the stress in the stringer. This formula 
is probably always conservative in the range in 
question. 

COMPARISON BETWEEN TEST AND CALCULATED RESULTS 

In order to check the validity of the method thus Ear 
developed, a test specimen was built to represent a 
structure corresponding to figure 1 (a). A sketch of the 
actual test specimen is shown in figure 5. Pin-end steel 
bars (not shown in the figure) spaced 3 inches apart 
were used to separate the edge stringers from the cen- 
tral stringer and to take up the transverse component 
of the diagonal-tension field that developed under load. 
In each bay between these bars, the strains in the 
stringers were measured with 2-inch Tuckerman strain 
gages on both sides of the specimen. This precaution 
proved necessary because the stresses on the two sides 
differed so much at some stations that readings on only 
one side would have been almost useless. 

The load was increased from zero to the maximum 
of 4,800 pounds in five steps. With a very few minor 
exceptions, the points for any one gage fell on straight 
lines. For each station, the results obtained on the 
front and the back of the specimen were averaged and 
the average values are plotted in figure 6. 



•5 6,000 




Figure 6. 
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Comparisons between calculated and experimental results for tension test panel. 
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Figure 1.— Comparison between calculated and experimental results for compression test panel. (Data from reference 2.) 
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The calculations were made for the two different 
assumptions of the shear stiffness indicated on the fig- 
ure. The second assumption of G e =%G in the top part 
was based on the experimentally observed fad that one 
well-developed diagonal-tension fold showed in the top 
of the panel on each side, in agreement with the cal- 
culation showing thai al the maximum load the shear 
stress in this region was about six times the buckling 
stress. 

The second assumption gives perfect agreement be- 
tween calculated and test results for the stress in the 
central stringer. The agreement is not quite so good 
on the edge stringer, the discrepancy occurring chiefly 
at the root. Several explanations of the discrepancy 
may be offered. An error of several percent may be 
caused by an error in the value of E assumed to convert 
strain readings to stress readings. The simple theory 
used may break down to some extent near the root and, 
finally, jig deflection may cause errors. The steel 
t riangle used on t he lower end is not a rigid foundation, 
and a slight elastic deformation of this steel triangle 
under the edge stringers would relieve the edge stringers 
of some load and throw it into the sheet and possibly 
into the central stringer. A deformation of about 
0.0003 inch would be sufficient to make the calculated 
st ringer stresses equal at the jig end. Undoubtedly the 
assumptions of effective areas, effective shear stiffness, 
and jig deflection could be varied within their possible 
limits to give a much better agreement with the experi- 
mental points. 

A similar analysis was made for the panel tested in 
compression as described in reference 2. The results 
are shown in figure 7. It will be noted that fair agree- 
ment with the experimental points is obtained by assum- 
ing that the effective shear stillness is only 0.2 the shear 
modulus, in marked contrast to the tension panel. The 
curves calculated with G e =G are also given to show 
the extent to which possible variations in G e affect the 
stringer stresses. 

BEAMS WITH ONE LONGITUDINAL 

BEAM OF CONSTANT DEPTH 

The simplest case of a beam subjected to shear defor- 
mation of the llange is shown in figure 8. For simplicity 
of the sketch the llange material on the side not under 
consideration i- assumed to be concentrated at the shear 
web. This assumption does not influence the analysis 
when the cover is flat. 

For convenience of discussion, the material concen- 
trated at the to]) of the shear web will he referred to as 
the "flange" throughout this paper, while the stringer 
attached to the cover sheet will be referred to as the 
"longitudinal.'' 

It is again assumed that the longitudinal is cut along 
the line of symmetry (fig. 8 (b)). The force acting on 



this halved longitudinal is denoted by F L , the force on 
the (tension) flange by F F . The shear force in the web 




(b) 



Figure 8.— Beam with flat cover and one longitudinal. 

is denoted by S w ) the shear force in the cover sheet, 
by S c . 

The governing equations are 

dx 

(I F F =Sw~j l — dS c (3 a) 

-dF L =dS c (3b) 
G 

dr= —ffifaF— *l) dx (3c) 

with the auxiliary equations 

<rp s =-j^; <?l~^; Sw=P; dS c =rtdx 

The solution of the resulting differential equation is 
given in appendix B, Case 3 (a). 

COMPARISON BETWEEN TEST AND CALCULATED RESULTS 

The test panel that had been used in the previously 
described tension test was slightly modified and 
attached to two duralumin I-beams to form an open 
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(b) Open side. 
Figure 9.— View of test beam, showing strain gages. 
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Figure 10. — Set-up for testing beams. 



1931 — 37 2 



s 



REPORT No. 608— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



box beam. Figure 9 shows photographs of the beam 
with the strain gages in place for a test run; figure 10 
shows the test set-up. The cross section of this beam 
is shown in figure 1 1 . 

It should be noted that the cover sheet and the longi- 
tudinal were not attached to the bulkheads except at 
the root. The flange material of the I-beams (includ- 
ing the cover strips riveted to them and the sheet 
material effective in tension) was replaced, for the pur- 
pose of analysis, by equivalent concentrated llanges 
with a centroidal distance of 2. SO inches (effective depth 
h of beam, fig. 8 (a)). The calculated stresses are 
therefore valid for the flange centroids. For compari- 
son with the measured stresses, the calculated flange 
stresses were corrected to the outside fiber stresses 



BEAM OF VARIABLE DEPTH 
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under the assumption that plane cross sections remain 
plane for the I-beams with cover strips. 

Figure 12 shows the experimental points, the curves 
calculated for three different assumptions of the shear 
stiffness, and the stresses calculated by the ordinary 
bending theory. It can be seen that the experimental 
points group fairly well about the curve for G e =% G, 
particularly when this curve is corrected for an esti- 
mated jig deflection by the formula in appendix B, 
case 2. Close to the root, however, discrepancies are 
again observed as in the case of the tension panel. 
The high flange stress at the station nearest the root 
may perhaps be explained by nonlinear stress distri- 
bution in the I-beams caused by the method of attaching 
them to the jig, which was not designed for this test. 
The reduction in shear stiffness of the sheet as compared 
with the stiffness developed by the same sheet in the 
tension panel can be ascribed to numerous initial 
buckles present in the beam but not in the tension 
panel. 

Inspection of figure 12 shows that very large varia- 
tions of shear stiffness have only a relatively small 
influence on the bending stresses. This result is due 
to the fact that, even when the shear stiffness increases 
to infinity, the bending stresses never exceed a finite 
limiting value. In many actual structures, the shear 
stiffness provided is sufficiently large to permit the 
limiting stress to be approached within a few percent. 
Practically speaking, this [act means that the shear 
stiffness need not be very accurately known to obtain 
the necessary accuracy in the bending stresses. 



In a beam with variable depth, the only change in the 
equations is introduced by the fact that the vertical 
components of the flange forces balance part of the 
applied shear, so that the shear in the web now becomes 



S w =S a -j (tan 0+tan y) 



(5) 



where p and y are the angles of inclination of the 
tension flange and of the compression flange. 

The analytical solution for a special case of a beam 
with variable depth is given in appendix B as Case 3 (b). 
4,000 z 



3,600 




0 6/3/6 34 30 

Distance from root, in. 

Figure 12.— Comparison between calculated and experimental results for test beam. 



CONSTANT-STRESS SOLUTION FOR BEAMS WITH ONE 
LONGITUDINAL 

The analytical solutions presented thus far, together 
with the trial-and-error method, are reasonably ade- 
quate for dealing with beams having one longitudinal. 
There appears to be but slight possibility, however, of 
extending these solutions to the practical cases of beams 
with a number of longitudinals. An approximate 
method will now be developed that can be extended to 
such beams. The method will first be developed for a 
beam with a single longitudinal because comparisons can 
be made with the exact solution to gain some idea of the 
reliability of the approximate method. 
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The approximate method is based on the following 
reasoning. It is the aim of the designer to dimension 
the structure so that the stress in it is uniform for the 
given loading. For several reasons this ideal is never 
reached, but there is usually an effort made to taper the 
dimensions so as to approach the dimensions of the ideal 
design. Now the solution for constant stress along the 
span can be very easily obtained. It is possible, there- 
fore, to consider the actual condition as a super- 
position upon the ideal case, which can be calculated 
exactly, of some additional disturbing cases or "faults." 
These faults can be calculated only approximately, but 
if they are of minor importance compared with the ideal 
case, the resulting error of the total solution will be 
small. 

The detailed development of the method is as follows: 
The fundamental equation 



G 



(6) 



can be integrated once, if o> and a L are constant as 
assumed, to give 

_ (o>— <r L ) C z n j (°>~ <ri)xGi 



Kb 



Eb 



(7) 



where G x is the shear stiffness averaged over the 
distance x = 0 to x = x, and the x origin is taken at the 
root. Integrated again to give the total shear force in 
the cover sheet 

S c = f 0 L r z tdx=K 1 (<T F -a L ) (8) 
For example, if G e and t are constant along the span, 

G e tU 



2Kb 



Equation (8) furnishes one relation between a F and a L . 
One more relation is needed to complete the solution. 
There are infinitely many conditions from which to 
choose this relation. At any station along the span, the 
internal bending moment should equal the external 
bending moment. The root section has been chosen 
because in a number of trials it always proved, by far, to 
be the best choice. Equating the internal and external 
moment (applied at the root) gives the relation 

(<t f A Fq + <y L A Lo )h 0 =M ao (9) 
Now remembering that 

Sc=ctlA Lq 

equations (8) and (9) can be solved for the bending 
stresses 



MpK, 

^UAf^Alo+KM^+A^)] 

M 0 (A Lo +KO 



(10a) 



(10b) 



* F holA^+KMFo+A^)] 
Substituting equations (10a) and (10b) into equation 
(7) gives 



■> <>' M 



H> + K i+ ^:)] 



(KID 



Equations (10a), (10b), and (10c) constitute the "pure 
constant-stress solution" for a beam with a single 
longil U( I inn 1 . 

The internal bending moment at any station alonu- 
the span can now be calculated 

M int = (a F A F + <r L A L )h 

and, in general, this internal moment will not be equal 
to the applied moment M a . This difference constitutes 
the first fault of the constant-stress solution and will 
be called the "moment fault." 

In order to remove this fault, additional (corrective) 
bending moments must be added, which arc at any 
station 

M' = M a -M int 

the prime denoting corrective moments. The stresses 
caused by these corrective moments must be computed 
and added to the stresses of the pure constant-stress 
solution. 

The method of computing the stresses caused by the 
corrective moments will be approximate and arbitrary 
as thus far no exact solutions of this problem have been 
found. The following method was chosen because the 
underlying assumption is the most obvious one and 
because the method is very convenient, eliminating the 
necessity of computing the internal moments, the cor- 
rective moments, and the corrective stresses separately. 

From equations (10a) and (10b) it follows that the 
ratio 



(ii) 



The assumption is now made that I his ratio remains con- 
stant (r=r 0 ) along the span and that it holds not only 
for the stresses caused by the "ideal" moments but also 
for the stresses caused by the corrective moments. 
Under this assumption, the direct stresses at any station 
are given by 

M a 



a F - 



(112a.) 



(12b) 



From these stresses the shear stresses axe obtained by 
using the fundamental relation (2) and integrating from 
the root toward the tip 

T=J* ( ^ h (<T F -<T L )dx (12c) 

The moment fault has now been removed ; that is, the 
internal moments equal the applied moments when the 
stresses as given by equations (12a) and (12b) exist in 
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the flange and in the longitudinal. But equation (12c) 
follows directly from equations I 12a ) and (121)) and the 
stresses given by (12a) and (12c) will not, in general, 
fulfil] the fundamental equation (3a) of equilibrium of 
the flange element. Equation (3a) requires that, for 
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Figure 13— Correction factor, G 
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equilibrium of the flange element, 
shear force in the cover should be 

AS cs = ^Ax-AF F 



the inclement of 



(13) 



where the additional subscript S denotes the increment 
required for static equilibrium. The incremenl of shear 
force actually developed is 

AS CE =rtAx (14) 

where the subcript E refers to the fact that this incre- 
ment is provided by the elastic deformations of the 
flange and the Longitudinal. Failure of the shear-force 
Increments given by equations (13) and (14) to be 
identical constitutes the second fault of the constant- 
stress solution, the so-called "shear fault." 

Static equilibrium for the flange elements would be 
restored if corrective shear-force increments were in- 
troduced equal to the differences of these two sets of 
shear-force increments 



AS c ' = AScs-*Sc 



(15) 



where the prime again denotes a correction. The cor- 
rect he shear force S c ' at any station is obtained by 
integrating from the tip to the desired station, the force 
being zero at the tip. The corrections to be added to 
the stresses would then be given by 



/ Sc 
- 1 / 



Sc 



AS c f 
' tAx 



(Care must be taken in determining the signs of the 
corrective stresses. The safest method is to compare 
their direction with the direction of the stresses given 
by the pure constant-stress solution.) 

Introducing these corrective stresses would restore 
static equilibrium but would again upset the basic 
elastic relation given by equation (6). A compromise 
must therefore he made by using only a fraction C x 
of the correction 



0> — Kj\ a 
- ' I- 



<tl — — C-i-j- 



AN,' 
IX, 



(16) 



These stress corrections are added to the stresses 
obtained from equations (12a), (12b), and (12c) to 
obtain the final corrected stresses o> . a L , and 

T corr- 

Values of C x may be established by comparing a 
number of exact solutions with the corresponding 
constant-stress solutions; an averaged curve is shown 
in figure 13. 

In order to gain some idea of the range of applicability 
of the constant-stress solution, a series of related beams 
was calculated. The characteristics of three of these 
beams are given in table II. The first set of calculations 
was made by using the analytical solutions given in 
appendix B for beam A and by using the trial-and- 
error method for beams B and C. The second set of 
calculations was made by using the constant-stress 
solution as described. The results of the calculations 
are shown in figures 14 to 16. 

For beam B, the stresses given by the pure constant- 
stress solution are also shown. Beam B is a constant- 
stress beam when analyzed by the ordinary bending 
theory and has zero moment fault. The complete 
analysis for this beam is given as an example in 
appendix C. 

It is to he expected that, in general, there will he 
smaller differences between the constant-stress solu- 
tion and the exact solution for beams with small moment 
fault than for beams with large moment fault. This 
expectation is borne out by the results. Beam 
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Figure 14.— Stresses in beam A. 



Tip 



which comes close to the ideal case, shows smaller 
differences than beam A, which is further from the 
ideal case because the areas A v and A L are constant 
along the span. Beam (\ which corresponds to an 
actual case, as far as variation of A Fi A L) /, and h along 
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the span is concerned, shows also good agreement for 
the bending stresses. The agreement is not quite so 
good for the shear stresses. 

Considering all the factors involved, it seems safe to 
assume that the constant-stress solution will give 
satisfactory results in practical cases for the maximum 
stresses, provided that the correction introduced by the 
shear fault is not larger than about 20 percent of the 
stress given by the pure constant-stress solution. 

BEAMS WITH MANY LONGITUDINALS 

YOUNGER'S SOLUTION 

Actual wing structures are built as box beams with 
many longitudinals, and the depth of the beam as well 
as all cross-sectional areas varies along the span. 

The first attempt at obtaining a solution for a multi- 
stringer beam was made by Younger (reference 4). He 
considered the limiting case of infinitely many longi- 
tudinals (i. e., a plate cover as shown in fig. 17) and 
assumed the box to be of constant section; for the dis- 
tribution of the bending moments he assumed a cosine 
Law. 

Younger's solution and its extension to arbitrary 
moment curves are given in appendix B. It should be 
noted that this solution does not fulfill the equation of 
equilibrium for the flange element (the differential 
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Figure 15.— Stresses In beam B. 



equation does not hold along the flange) so that a shear- 
fault correction is necessary, as discussed in connection 
with the constant-stress solution for the beam with a 
single longitudinal. 



The usefulness of Younger's solution is so limited by 
the assumption of constant cross section along the 
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Figure 16.— Stresses in beam C. 



span that a more general method appeared desirable. 
The constant-stress solution was developed to fill this 
need of practical stress analysis. 

The principles of the constant-stress solution have 
been discussed in detail for beams with a single longi- 
tudinal. The extension of the solution to beams with 
many longitudinals is given in appendix B. The 
practical procedure of applying it is essentially identical 
with the procedure outlined for beams with a single 
longitudinal. The constant K 2 is computed and used 
to compute the constant K 3 for the root section, using 
equation (B-27). The stresses at a number of stations 
along the span are then obtained by the formula 



M cosh K z y 

h(A F cosh K 3 b+-j^ sinh K z b^ 



(17) 



where y varies from y—0 for the center line of the beam 
to y =6 for the flange. The shear stress in the cover 
sheet next to the flange is obtained by integrating from 
the root outward the expression 



V dx/ F 



a F K s tanh K 3 b 



(18) 



where o> is obtained from equation (17) by setting 
y=b. Equation (18) is obtained from equations (B-20) 
and (B-25). 
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The increments of corrective shear force are obtained 
by using equations (13), (14), and (15). After the 
integration of ( If)) in from the tip to obtain the correc- 
tive shear force S c ' y the correction to the flange stress 
is calculated by the first expression of (16); the correc- 




Figure 17.— Notation used for beams with orthotopic cover plates. 

tion to the shear stress is calculated by the last ex- 
pression of (16). 

The calculation of the correction to the stress a L is 
somewhat more complicated because it varies along 
the chord. The total force on all longitudinals, using 
equation (17), is given by 

Fl=J* *T dy =K£ ffCL sinh Kzh (19) 

where a C L denotes the stress at the center line of the 
beam obtained from equation (17) by setting y=0. 
In accordance with (16), only a part of the corrective 
shear force is applied so that the corrected total force 
on the longitudinals is 



I r L corr =F L —C l Sc r 



rjn. 



Assume now that the corrected stresses in the longi- 
tudinals are distributed chordwise according to the law 



(21) 



(22) 



<r C orr=(rcL corr cosh Yy 
The unknown Y can be found from the equation 

tanh Yb _ F LcoTT 
Yb ~A L <r Fcorr 

which is based on the premise that 

^ L corr ' corr 



for y= b. After Y has been found, the corrected stress 
at the center line is found from 



sech Yb 



and equation (21) can then be used to calculate the 
stresses at intermediate values of y. The right-hand 
side of equation (22) is the ratio of the average stress 



in the longitudinals to the stress in the flange. In 
general, this ratio will be less than unity; however, 
figure 16 shows that for a beam with a single longi- 
tudinal the stress in the longitudinal may be larger 
than the stress in the flange over a part of the span, and 
similarly the right-hand side of equation (22) some- 
times may exceed unity. In such a case, equations (21) 
and (22) may be replaced by 



Vcorr=o-cL corr (2 - cosh Yy) 

/ sinh Yb \ 
V Yb J _J_£co 
(2 -cosh Yb) A L a Fi 



corr 
corr 



(21a 



(22a 



After Y has been found, the corrected stress at the 
center line is found from 



(2 -cosh Yb) 



and equation (21a) can then be used to calculate the 
stresses at intermediate values of y. 

The solution of equations (22) and (22a) can be 
effected by inspection of tables. Kor practical pur- 
poses it should be sufficient to use the curve given on 
figure 18. 

As examples, beams A and B were analyzed under 
the assumption that longitudinals with the total cross- 
sectional area A L are distributed uniformly along the 
chord. The results are shown in figures 19 and 20. 
It will be seen that the stress at the center line of the 
beam is very low. If all longitudinals are of the same 
cross section, they must be designed to the stress in the 
first longitudinal adjacent to the flange. Consequently, 



r \ 




1 1 






\ 










i 


.0 








'-corr 








1 
1 














8 


s 












1 
l 


























1 
I 












.5 
b 

.4 


6 














— 1- 
1 

— 1— 






i 




—1 


Yb 

-4— 




\ c 


Ise scalt 
it left 


» 




1 
l 


























1 
1 

— U 








se s 
t ri 


co/e 
ght 
















1 
1 












.1 


2— 














I- 












■0 — 


























-0- 



Ratio 



Ar(T f 



+ U F CO rr 



Figure 18.— Graph for auxiliary parameter Yb. 

the longitudinals near the center line are very in- 
effectively used. In this connection, attention might 
be called to the fact that the longitudinals need not be 
of the same cross-sectional area along the chord. The 
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assumption of A L being uniformly distributed may be 
fulfilled, for instance, by using longitudinals of large 
cross-sectional area but widely spaced near the flange 
and longitudinals of small cross-sectional area but 
closely spaced near the center line. Although such an 
arrangement would not increase the over-all structural 
efficiency, it might under certain conditions offer 
man ufacturing advantages. 

MECHANICAL ANALYZER 

The constant-stress solution is always approximate. 
When the moment and shear corrections are large, 
doubts may arise as to whether the solution is suffi- 
ciently accurate. It might be advantageous to con- 
struct a mechanical analyzer to deal with such cases. 
One possibility for such an analyzer would be actually 
to build units representing the mechanical model 
sketched in figure 1 (b). The springs might be canti- 
lever springs, so that their stiffnesses could be varied by 
changing their lengths. Each unit would represent 
one bay of the trial-and-error method of solution and 
would have one spring to represent the stringer stiffness 
and one spring to represent the shear stiffness of the 
sheet attached to one side of the stiff ener. 

The chief difficulty in the design of such an analyzer 
would probably be in reducing the friction between the 
units and the guides necessary to aline them. A fairly 
large number of units would be necessary to represent 
a wing cover, which would mean a fairly expensive 
instrument. This disadvantage is counterbalanced by 
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Figure 19.— Stresses in beam A with Al uniformly distributed along chord. 



the possibility that the instrument would offer in a 
comparatively short time quite an exact analysis, 
including the effects of bulkheads and of yielding 
supports. The main errors in this solution would be 
those caused by the finite length of bays. 



CONCLUSION 

The art of stress-analyzing shell structures is of recent 
origin, and any methods of analysis proposed must go 
through a process of trial and development. 

Development of the method of shear-deformation 
analysis is desirable in several directions; e. g., exact 
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Figure 20.— Stresses in beam B with A L uniformly distributed along chord. 



solutions should be found to replace the constant-stress 
solution and methods should be devised to calculate 
t he influence of bulkheads. 

Rough approximate calculations on bulkhead effect 
can be made by assuming that all the longitudinals 
are relocated at the center line of the beam. For 
beams with a single longitudinal, the effect of bulk- 
heads can be calculated. A series of systematic com- 
parisons between the extended solution of Younger and 
Case 3 (a) of appendix \\ indicates that for a certain 
range the single-longitudinal assiimpt ion may yield 
acceptable approximations when used in conjunction 
with suitable correction factors. The comparisons 
are not given, however, because they might be mis- 
leading in view of the shear fault of Younger's solution. 
Calculations made thus far indicate that in practical 
cases the effect of the bulkheads is very small. 

It should be emphasized that analyzing shell struc- 
tures is an art rather than a science. The arithmetic 
of analyzing highly redundant structures can be re- 
duced to manageable proportions only by making 
assumptions that will be valid only within a certain 
range. This fact leads to the unfortunate, but inevi- 
table, conclusion that the analysis of such structures 
cannot be made entirely by handbook and formula but 
must be guided by engineering judgment. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 3, 1937. 



APPENDIX A 



LIST OF SYMBOLS 

A, cross-sectional area (sq. in.). 

E, Young's modulus (lb. per sq. in.). 

F, internal force (lb.). 

G, shear modulus (lb. per sq. in.). 
K, constant. 

L, length of panel or beam (in.). 
.1/. bending moment (in. -lb.). 
P, external load (lb.). 
S, shear force (lb.). 

/>, spacing of stringers (in.). (See figs. 3 and 4.) 

b, half width of beam (in.). (See fig. 8.) 

c, camber of cover (in.), 
depth of beam (in.). 

t, thickness of cover sheet (in.). 

ti, displacement of point (in.). (See fig. 4.) 

w, running load (lb. per in.). 

y, shear strain. 

<Tj direct (normal) stress (lb. per sq. in.), 
r, shear stress (lb. per sq. in.). 
14 



Subscripts have the following significance: 

A, loaded stringer A shown in figures 1, 2, 21, and 22. 

B, unloaded stringer B shown in figures 1, 2, 21, and 22. 

C, cover sheet. 

/•'. flange of beam. 

L, longitudinal of beam. 

U'. shear web. 

a, applied shears and bending moments. 

r, effect ive. 

0, root section. 

c, compression. 

/, tension. 

int, internal. 

corr, corrected. 

S, static equilibrium. 

/•'. elasl ic equilibrium. 

CL, center line. 



APPENDIX B 



SOLUTIONS OF DIFFERENTIAL EQUATIONS FOR SYMMETRICAL STRUCTURES OF CONSTANT 

CROSS SECTION 



SIGN CONVENTIONS 

Forces and stresses in stringers are positive when 
tensile. Shear forces and stresses in the sheet are posi- 
bive when caused by positive stresses a n< I strains in the 
loaded stringer A in the case of axially loaded panels 
or in the flange F in the case of beams. 

CASE 1— THREE-STRINGER PANEL ON RIGID FOUNDATION WITH 
AXIAL LOAD 

The two possible cases shown in figures 21(a) and 
21(b) can be mathematically treated by taking one- 
half the panel, as shown in figure 21(c), which also 

////////////////////// ///////////////////// ///////////////// 



-b 



(a) (b) (c) 

Figure 21.— Axially loaded panels. 

gives the notation to be used. The derivation of the 
l'u iu la mental equations is given in the main body of 
this paper. Slightly modified for the purpose of deriv- 
ing the basic differential equation, these equations are 



0"A 



--T- and a B =— -r- 



Ge_ 
Kb 



(B-l) 



(B-2) 



w here the primes denote dhl'erent iation with respect to J*. 

Differentiating equation (B-2) again and substituting 
into the result from equation (B-l), 

The boundary conditions are 
at x=0 , r=0 



a I 



The result i- 



P 



P_ _Ge_ sinh Kx 
~~A A EbKcosh KL 



(B-4) 



As 

A a 



30 



cosh Kx \ 
cosh KL/ 



w here 



* 2 =§(i+i) 



I u reference 2 I he formula 



, 2Pr cosh px— tanh pL sinh px , ~| 



u here 



p .-25 G e t 
v ~~ 2 A T hE 



is given for the special case where the area of the edge 
stid'ener is twice the area of the cen l ra 1 -tiilener. Tak- 
ing account of the (inferences in notation and coordinate 
systems used, this result agrees with the general formula 
given under (B-5). 

It should be noted that the final formulas (B-5) be- 
come invalid when either / or G e approaches zero be- 
cause in these cases the equation (B-3) becomes invalid. 
The solution for such cases is obtained by using the 
fundamental equations (B-l) and (B-2) directly. 

An analogous procedure must be used for Cases 2 
and 3. 

CASE 2— THREE-STRINGER PANEL STRAINED BY MOTION OF 
SUPPORTS 

The differential equation for the case of figure 22 is 



(a) (b) fc) 

Figure 22.— Panels strained by motion of supports. 

the same as for Case 1. The boundary conditions are 
now : 

at x=0 , a A = 0 and a B = 0 ' 



i\ t 



x—L } t= j G e =T 0 



(B-7) 



The result is 

cosh Kx 
r ~ ro cosh KL 

t sinh Kx 

t sinh Kx 
° B ~~ t °KAb cosh KL 

where K has the same meaning as in (B 6). 



(B-8) 
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CASE 3— CANTILEVER BEAM WITH ONE STRINGER 

(a) Uniform depth, concentrated load at tip. 
(6) Depth decreasing lineally to zero, uniformly dis- 
tributed load. 

Figure 8 shows the notation used for both cases. 
(Note that the x origin is at the tip.) The funda- 
mental equations are for Case 3 (a) 



<r/A F 



h 



-rt 



a L 'A L = Tt 



f G c , 

7 — M {aF ~ 



l>> 1) 



which gives the differential equation 

PG e 



ApEbh 



=0 



The boundary conditions are 



at 
at 

The result is 



x=0, a F =Q, and o L - 
x = L, t = 0 



(B-10) 



(B-ll) 



thl 



cosh Kx \ 
' cosh KLJ 



sinh Kx ^ 
K cosh KLj 



(B-12) 



CL= h(A L +A r ){ x ~ 

w r here K has again the same meaning as in (B-6) with 
A F and A L substituted for A A and A B . 

In Case 3 (6), wL/2 is substituted for P; h in this 
case is the depth at the root. 




2* 



Figure 23.— Cantilever beam with concentrated load not. at tip. 

The case of a beam loaded by a concentrated load 
not at the tip is a simple problem in indeterminate 
structures. The beam is cut just outboard of tbe load 
(fig. 23) and the stresses in the cantilever pari are cal- 
culated (Case 3 (a)). From these stressss, the distortion 
of the beam section at the cut; L e., the relative dis- 



placement of the tips of the flange F and the longi- 
tudinal L, can be calculated. A system of forces X 
is then applied to equalize the distortion of the can- 
tilever tip and of the inboard end of the "overhang," 
utilizing the formulas of Case 2. 

CASE 4— CANTILEVER BEAM WITH ORTHOTROPIC COVER PLATE 

Younger's solution for a beam of constant section.— 

The beam and the coordinate system used are shown in 
figure 17. It should be noted that the x direction is 
opposite to that used in Cases 3 and / h 

Under the assumptions that the transverse stresses 
and strains are negligible (Poisson's ratio equal fco zero), 
and that G e is independent of A\ the differential equa- 
tion of the cover is 



(B-13) 



where u is the displacement of any point on t he cover in 
the ./• direction. 

The boundary conditions are 

bu 

at x=0, u=0 and 5— =0 

by 



x=L, ^— = 0 

Or 



(B-14) 



y- 



by 



This equation was established by Younger (reference 
4, pp. 36-47). For the solution he assumed thai the 
external bending moment (on the whole beam) is 
given by 

M=M 0 cos (B-15) 
and obtained for the longitudinal stress in the cover 



M 0 cosh cos 



2/(1 A F cosh 



7T& 



IK I. 

and for the shear stress 



2KL 

7T 



sinh /kl) 



sinh 



Try . tx 
2ivZ Sin 2l 



' hE nTrA , x6 . AKHL . , irb 
2KA F cosh 2^Jj ~ — smn 



1KL 



(B-16) 



where K is defined by 



JS2 

A ~E 



(B-17) 



Extension of Younger's solution. -Younger's solution 
can be somewhat extended. The external bending 
moment can be represented by a superposition of 
several terms: 



\ r w 7TCC . -m . ''Wx . . i 07r.r 
M=M X cos 2X+M3 cos 2^+M 5 cos 



+ 



mirx 



. f in n .1 

+ M,„ COS 7f£ 



<B 18) 



where t he m's are odd integers. 
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The values M\ . . M m are chosen so that the sum of 
the terms equals the given external bending moment at 
m points other than the tip, where it is assumed that 
M=0. In order to make comparisons with Case 3, the 
bending moment caused by a tip load was expressed by 

M=PL ^0.821 cos 11+0.101 cos |g 



+ 0.045 cos ^+0.033 cos 



7tvx\ 
1L I 



(B-19) 



The stresses corresponding to the rath term are given 



by 



m wx 



M m cosh m cos w 



2h( K A F 



cosh 



mwb , A L 2KL . , m7r&\ 
2KL+T -^ smh 2KL) 



Tm ~ h E 



. . in 7t ii . mirx 
3wh 2KL am ^L 



( 2 



2KA F cosh 



trncb , ArAK 2 L 



1KL 



himv 



smh 2KLj 



(B- 
16a) 



The assumptions of Poisson's ratio being zero and 0 
being independent of A' are, strictly speaking, incom- 
patible. The physical picture conforming to these 
assumptions is not a plate but a system of stringers 
carrying only longitudinal stresses tied together by a 
sheet carrying only shear stresses. This picture is 
realized very nearly in practice by a skin-stringer cover, 
the only difference being that the total cross-sectional 
area of the stringers is not necessarily equal to the area 
of the sheet, as in the case of the plain cover sheet. All 
the equations written for the plain cover sheet apply, 
therefore, to the skin-stringer cover if only (B-17) is 
replaced by 

K 2 =rQ (B-17a) 

where R is the ratio of sheet, area to area of longi- 
tudinals. 

Constant-stress solution. — The coordinate system is 
thai sihown in figure 17. ruder the assumption that 
o-= constant for each longitudinal, the fundamental 
relation 

H£ (B - 20) 

can be integrated once to give 

Ao- p r , xAa ^ 



(B 21 



w here G x is the shear stiffness averaged over the distance 
x=0 to x=x. Integrating again 

Sc=f o L Tjdx=f o L ^G x t</, 



In any given case this integration can be performed 
and the result is 



S C =K 2 



where 
Now 

(see fig. 24) or 



-\'/ 



(H 22) 



K 2 = jgGjdx 



dy " b 



(B-23) 




A F £ 

Figure 24 —Free-body diagram of cover plate. 



I differentiating (B-22) and equating to (B-23) 



d 2 a A L 



T hK. 



(B 24) 



assuming that K 2 is independent of y. 
The boundary conditions are 



da 



(1) at y=0, t=0 for any x. Therefore 

(2) at any desired reference station R, the internal 
moment equals the external moment M R . 

The solution is 



M R cosh K%y 



hJ^<A F cosh KJ) + ^sm\\Kzb^ 

~Q 

T =<7^K z x tanh K 3 y 
where K 3 is defined by 

rs 2 ^£ 



(H 2.r) 



(B-26) 



(B-27) 



-gG x tdx 



It may be noted that if G e and t are not varied along 
the span, the constant K 3 is identical will) the corre- 
sponding constant of Younger's solution except for a 
10 percent difference in the numerical factor, namely, 
-y/2 against w/2. 



APPENDIX C 



ANALYSIS OF BEAM B 

The dimensions and the loading of the beam are 
shown in table 1 1 . 

ORDINARY BENDING THEORY 



M 



2,800,000 



<j F -o L -G- h ( ^ + A ^ 24(1.875+1 .875) " 
31,100 lb. per sq. in. 

CONSTANT-STRESS SOLUTION 

Since G e is assumed constant along the span, G x =G e 
and, from equation (7), 

From equation (8) 



(c 

- tLO Jo \ 

,0 .2X0.040 f 280 V, x \ , 

— 24— jo \ 1 -mr t 

=4.35 (o>— (t l ) 
^1=4.35 

From equation 1 1 1 la 

= 2,800,000X4.35 

° L 24[1.875X1.875+4.35(1.875 + 1.875)] 
=25,550 lb. per sq. in. 

From equation (10b) 

2,800,000(1.875+4.35) 



24[1.875X 1.875+4.35(1.875 + 1.875)] 
= 3(i, 500 lb. per sq. in. 
18 



Substituting in equation (7) for the shear stress at 
the tip 

Tmax= (36,500-25,550) ^ ' =25,560 lb, per sq. in. 

The calculation of the shear correction is shown in 
table III. 

TRIAL-AND-ERROR SOLUTION 



Take Az=40 in. 
S w Ax wxLAx_ wL 



71.4X280X40 1AA(ynll 
Aj= kt^-fta =16,6/0 lb. 



h ' 2hoX ~2h^~ 2X24 
AF F =16,670-AS C 
GAx, x 0.2X40 

~Eb~ ^ = — 24 — ^ F ~~ * L ' =0 - 333 c u 
A typical cycle of the calculation is shown in tabic I V 
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Aa=0. 403 
Ab=0. 220 
*=0. 016 
6=4.60 



L=36 
Ax= 6 

Fa=2, 400-2ASc 

Fa Fa 
aA ~AA 0.403 



Fb=2ASc 

Fb Fb 
° B Ab 0.220 
G e /E=(}.4 



ASc- 



= %ft £ (<M-<r«) =0.522 (<r A -o H ) 
= 2 At 

-t(Ax=0. 096 r 



Station 


By trial-and-error method 


By formula 1 


ASc 
(lb.) 


Fa 
(lb.) 


OA 

(lb./sq. in.) 


F B 
Ob.) 


OB 

(lb./sq. in.) 


a a— OB 

(lb./sq. in.) 


At 

(lb./sq. in.) 


T 

(lb./sq. in.) 


ASc 
(lb.) 


OA 

(lb./sq. in.) 


OB 

(lb./sq. in.) 


T 

(lb./sq. in.) 


0 




2,400 


5, 960 


0 


0 










5,960 


0 


5,230 




376 






3, 887 


373 


2,024 


5, 020 


376 


1,708 


3,312 


1,730 


5, 022 


1,717 


2, 885 


2 


210 


2, 157 


214 


1,814 


4,500 


586 


2,662 


1,838 


960 


4, 502 


2,670 


1,584 


3 


112 


1,197 


115 


i, 702 


4,224 


698 


3, 170 


1,054 


550 


4,220 


3, 186 


856 


4 


60 


647 


62 


1,642 


4,075 


758 


3, 444 


631 


329 


4, 070 


3, 461 


442 


5 


29 


318 


30 


1,613 


4, 005 


787 


3, 575 


430 


224 


3, 996 


3,595 


187 


6 


9 


94 


9 


1,604 


3,980 


796 


3,618 


362 


189 


3. 968 


3, 630 


0 











' Appendix B, Case 1. 



Table II.— CHARACTERISTICS OF BEAMS 

The beams are assumed to be half beams as shown in fig. 8 (a). 
All beams: 

h =24 in. at root. 6=24 in. 

/i=0attip. L=280in. 
G./E=0.2. 11=71.4 lb./in. 



Beam 


A F =A L 
(sq. in.) 


t 

(in.) 


c 


Root 


Tip 


Root 


Tip 


A 


1.875 


1.875 


0. 040 


0. 040 


0. 


B 


1.875 


0 


.040 


.000 


0. 


C. 


1.880 


.470 


.040 


.010 


0. 



Table HI.— CALCULATION OF SH KAR-F Al'LT CORRECTION FOR BEAM B 



Aj=40 

/V = o>vU = 36,500. If 
At 

AScs^Sw^-AF 



t= 25,560^ 

ASc*=t<A* 
ASc' = ASc5-AScif 



o F =0.5— 

^' = -0.5-^ 
_ =36,500+<r/ 



<r L =25,550+a// 

t =0.5—-— 
(Ax 

Teorr = T+T' 



Station 


i 

from 
root 
(in.) 


Ap =Al 
(sq. in.) 


h 

(in.) 


t 

(in.) 


O W -7- 
h 


F F 
0b.) 


AF 
(lb.) 


AScs 
(lb.) 


(lb./sq. 
in.) 


AScb 
(lb.) 


ASc' 
(lb.) 


Sc' 
(lb.) 


op' 
(lb./sq. 
in.) 


OL 

(lb./sq. 
in.) 


'r eorr 
(lb./sq. 
in.) 


(lb./sq. 
in.) 


(IbVsq. 
in.) 


Teorr 

(lb./sq. 
in.) 








































I. 


260 


0.134 
.268 
.402 
.536 
.669 
.804 
.937 
1.072 
1.205 
1.340 
1.473 
1.608 
1.740 
1.875 


1.71 


0. 00286 


18, 600 




9,800 


6, 860 


23, 740 


2,720 


4, 140 












18, 166 


41,840 


9, 800 


4, 140 


7, 720 


-7, 720 


44,220 


17, 840 


2 


220 


5. 14 


. 00857 


16, 660 


D.SOI) 


6,860 


20,080 


6, 880 


-20 


-30 


20, 050 


19, 600 


4, 120 


3,840 


-3,840 


40, 340 


21,720 


3 


180 


5. 86 


.0143 


16, 660 


9,800 


6, 860 


16, 440 


9, 400 


-2, 540 


-2, 220 


14, 220 


29, tOO 


1.580 


98G 


-980 


37, 180 


24, 580 


4 


140 


12.66 


.0200 


16,660 


9,800 


6,860 


12, 780 


10, 220 


-3, 360 


-2, 100 


10,680 


39, 200 


-1,780 


-830 


830 


35, 670 


26, 390 


5 


100 


15. 42 


.0257 


16, 660 


9, 800 


6, 860 


9, 130 


9, 400 


-2, 540 


-1,230 


7,900 


49,000 


-4, 320 


-1,610 


1,610 


34,890 


27, 170 


6 


60 


18. 87 


.0314 


16, 660 


9, 800 


6,860 


5, 480 


6, 880 


-20 


-10 


5, 470 


58, 800 


-4,340 


-1,350 


1,350 


35, 150 


26, 910 


7 


20 


22. 28 


.0371 


16,660 


9,800 


6, 860 


1,826 


2,710 


4, 150 


1,400 


3, 226 


68, 600 


-190 


-50 


50 


36,450 


25, 610 



























Table IV.— TRIAL-AND-ERROR SOLUTION FOR BEAM B 



Station 


x from 
root (in.) 


Ap=Al 
(sq. in.) 


tAx 


ASc 
(lb.) 


F L 
Ob.) 


ol 

(lb./sq. in.) 


AFp 

(lb.) 


F F 
(lb.) 


OP 

(lb./sq. in.) 


op—OL 

(lb./sq. iu.) 


At 

(lb./sq. in.) 


T 

(lb./sq. in.) 


ASc 
(lb.) 


0 


280 
260 
240 
220 
200 
180 
160 
J 40 
120 
100 
80 
60 
40 
20 
0 


0 






0 


0 




0 


0 


0 








1 


0. 1142 


3,730 


12, 940 




32, 896 


3, 755 


.268 


3, 730 


13, 920 


12, 940 


48,280 


34, 360 


11,453 


2 


. 3426 


7, 440 


9, 230 


21,443 


7,340 


.536 


11, 170 


21, 740 


22,170 


41,400 


19,660 


6, 553 


3 

4 


.571 


8,680 


7, 990 


14,890 


8, 500 


.804 


19, 850 


24,700 


30, 1G0 


37,500 


12,800 


4.267 


.799 


8, 670 


8,000 


10, 623 


8, 490 


1.072 


28,520 


26, 620 


38, 160 


3. r >, 0(H) 


8,980 


2, 993 


5 


1.028 


8,070 


8, 600 


7, 630 


7,840 


1.340 


36, 590 


27, 300 


46, 760 


34, 880 


7, 580 


2,527 


6 


1.255 


6, 360 


10,310 


5, 103 


6,400 


1.608 


42,950 


26, 700 


57, 070 


35, 470 


8, 770 


2, 923 


7 


1. 485 


3, 130 


13, 540 


2, 180 


3,240 


1.875 


46,080 


24,560 


70, 610 


37,640 


13,080 


4,360 
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V 

z 

Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designation 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal 


X 
Y 
Z 


X 
Y 
Z 


Rolling _ 


L 

M 

N 


y >Z 

z — >x 

X >Y 


Roll.. __ 


d 


u 

V 

w 


V 
Q 
r 


Lateral 

Normal 


Pitching 

Yawing 


Pitch 

Yaw 







Absolute coefficients of moment 
r _ L _ M 

Vl ~qbS Vm ~qcS 
(rolling) (pitching) 



D, 

P, 

PIT), 

V, 

V s , 

T, 

Q, 



Diameter 
Geometric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 

Thrust, absolute coefficient C T = 



Angle of set of control surface (relative to neutral 
n __jN_ position), 8. (Indicate surface by proper subscript.) 

Cn ~qbS 
(yawing) 

4. PROPELLER SYMBOLS 

P 

P, Power, absolute coefficient C P 
C 8y 



'pn'D 5 



Pn 2 



P n 2 D* 



Torque, absolute coefficient ^Q = ' p ^]j5 



Vy 

71, 



Speed-power coefficient = 

Efficiency 
Kevolutions per second, r.p.s. 

Effective helix angle = tan" {t^—} 



5. NUMERICAL RELATIONS 



1 hp. =76.04 kg-m/s=550 ft-lb./sec. 
1 metric horsepower =1.01 32 hp. 
1 m.p.h. = 0.4470 m.p.s. 
1 m.p.s. = 2.236l) m.p.h. 



1 lb. = 0.4536 kg. 

1 kg=2.2046 lb. 

1 mi. = l,609.35 m=5,280 ft. 

1 m=3.2808 ft. 



